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Abstract 

General scaling arguments, and the behavior of the thermal entropy density, are shown to lead 
to an infrared metric holographically representing a compressible state with hidden Fermi surfaces. 
This metric is characterized by a general dynamic critical exponent, z, and a specific hyperscaling 
violation exponent, 9. The same metric exhibits a logarithmic violation of the area law of entan- 
glement entropy, as shown recently by Ogawa et al. ( |arX iv:llll.lQ23|). We study the dependence 
of the entanglement entropy on the shape of the entangling region(s), on the total charge density, 
on temperature, and on the presence of additional visible Fermi surfaces of gauge-neutral fermions; 
for the latter computations, we realize the needed metric in an Einstein-Maxwell-dilaton theory. 
All our results support the proposal that the holographic theory describes a metallic state with 
hidden Fermi surfaces of fermions carrying gauge charges of deconfined gauge fields. 
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I. INTRODUCTION 



Much recent work [TH40] has focused on the holographic description of compressible metal- 
lic states of quantum matter. 

From a condensed matter perspective, a key issue in the zero temperature theory of any 
metallic state is the fate of the Luttinger relation [29j , equating the total charge density, Q to 
the volumes enclosed by the Fermi surfaces. An important feature of the early theories 

of holographic metals was that the Luttinger relation was badly violated: the holographic 
theory exhibits Fermi surfaces of gauge-neutral 'mesinos' (in the terminology of Ref. [39]), 
but they enclose volumes which are much smaller than the charge density, Q, of the boundary 
theory. (In the AdS/CFT correspondence, Q diverges with large numbers of colors, A^, while 
the mesino Fermi surfaces enclose volumes of order unity.) It was argued in Refs. [27[[29] that 
this deficit must be made up by hidden Fermi surfaces of 'fractionalized' fermions carrying 
gauge charges; we will refer to these gauge-charged fermions generically as 'quarks' (although 
they need not co-incide with the elementary fermions of any weak-coupling formulation of 
the boundary theory; in the context of supersymmetric gauge theories, we could also refer 
to them as 'gauginos', but we will not do so here). The holographic theory is only able 
to access gauge-invariant observables, and so direct signatures of the quark Fermi surfaces 
remain hidden to probes on the boundary. 

In subsequent holographic studies, states obeying the standard Luttinger relation were 
indeed found JEJ [171 [3ll [36] . However, they all had Fermi surfaces of gauge-neutral mesinos 
alone, and are all ultimately expected to be confining Fermi liquids (FLs) at low temper- 
atures. Our interest here is primarily on non-zero density states in which the gauge field 
is deconfined. Such states are expected to have Fermi surfaces of quarks alone, as in the 
'non- Fermi liquid' (NFL) states of Ref. [29j. Another deconfined state is the 'fractionalized 
Fermi liquid' (FL*) [ 4T[ 142] . which has co-existing quark and mesino Fermi surfaces [29j; 
holographic phases with charge fractionalization were considered in [39]. The key prob- 
lem of holographically detecting the hidden quark Fermi surfaces in these deconfined states 
remained open. 

In a recent paper, Ogawa, Takayanagi, and Ugajin [40j have proposed an elegant solution 
to this conundrum: compute the holographic entanglement entropy of the compressible state, 
and match it to that expected from the hidden quark Fermi surfaces. They considered a 
class of Einstein-Maxwell-dilaton (EMD) theories of compressible quantum states pTl [T3l 
[T4l [T9l [221 [231 [25] . and described conditions under which the area law of entanglement was 
logarithmically violated. The quark Fermi surfaces are expected |l3l - [i5] to display such a 
logarithmic violation of the area law, and so such EMD theories can be regarded as effective 
holographic descriptions of deconfined compressible states with quark Fermi surfaces, even 
though the EMD theories have no explicit fermionic degrees of freedom. Indeed, even in 
quantum models with only bosonic degrees of freedom, non-superfiuid compressible phases 
at non-zero Q are expected to have Fermi surfaces of fractionalized fermions [46] . 
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A. Violation of hyperscaling, entropy, and entanglement entropy 



Here we present a separate and simple general scaling argument for the infrared (IR) 
metric of holographic theories of non-Fermi liquid states with hidden Fermi surfaces in d 
spatial dimensions. Condensed matter physicists can find an introduction to the relevant 
concepts of the holographic method in Ref. [37]. We can write an arbitrary, holographic 
infrared (IR) metric in the scaling form 

*^ = ^(-;si;^ + '-"""'*^ + (1-1) 

where t is time, Xi {i = 1 . . . d) are the d spatial directions, and r is the emergent holographic 
direction. All numerical pre-factors (which are independent of t, x^, r) have been set equal 
to unity, and we have used reparametrization invariance in r to fix the co-efficient of dx'f to 
equal Such metrics can be realized for compressible phases at generic charge densities 

Q, and the numerical pref actors then depend upon Q in a manner which will be discussed 



in Section IL The boundary quantum theory is at r = (where the form (1.1) will not 
apply), and its low energy physics is captured by the IR limit, r ^ oc. This metric is 
characterized by two independent exponents, z and 0, which have been chosen so that scale 
transformations by a factor ( have the simple form 

t ^ Ct 

ds C^/Us. (1.2) 

This transformation makes it clear that z is the dynamic critical exponent. Relativistic 
conformal field theories have z = 1 and ^ = 0. Theories with 9 = but z ^ 1 are often 
referred to as "Lifshitz" theories j47l [48] . 



The novel feature of the metric (1.1) is that the proper distance ds of the emergent 
spacetime transforms non-trivially under scale transformations, with the exponent 9. In 
the usual AdS/CFT correspondence, the proper distance is invariant with = 0, rather 
than transforming covariantly under scale transformations. There is a natural connection 
between volume elements in the holographic space, and various entropic measures of the 
boundary theory, and so this suggests that a non-zero 9 will modify the scale transformation 
of the thermal entropy density, S. Indeed, using this reasoning we now show that 9 is 
the hyperscaling violation exponent [49j for the boundary theory; thus, non-invariance of the 
proper distance in the holographic theory implies violations of hyperscaling on the boundary. 
Hyperscaling is the property that the free energy scales by its naive dimension; specifically, 
in theories with hyperscaling has the temperature, T, dependence S ^ T^^^ [50l IST]. In 
the present holographic context, we can determine S by its proportionality to the area of 
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the horizon of a black brane which appears for T > 0; so for a horizon at r = r/^, we have 
from (1.1) that S ^ r^^. Also, from and (1.1) and (1.2) we determine that r scales as 



{d-e)/d, 



(1.3) 



and so we deduce that scales as t^^ ^^1^ . Noting that T is an inverse time, we can conclude 
that 

S ^T^^-^^l\ (1.4) 
establishing that Q is the hyperscaling violation exponent; an explicit computation of the 



entropy density obeying (1.4) will appear in Section II A, and analogous computations have 



appeared in other works [HI [I3l [Mj [25] . We are interested here in compressible states with 
fermionic excitations which are gapless on a (rf— l)-dimensional surface in momentum space 
(the Fermi surface) , and which disperse along the single dimension transverse to the surface 
with dynamic critical exponent z\ such states should have an entropy density 



S T"^/^, for general d\ 



:i.5) 



the entropy is that of the transverse excitations at each Fermi surface point, times the area 
of the Fermi surface. This is also the behavior found in recent gauge theories j46l [52l - [54] 
of non- Fermi liquid states in = 2, provided we identify z with the dynamic scaling of the 
Green's function of the gauge-charged fermions in the direction normal to the Fermi surface 
(in Ref. [53], the symbol z was used for the dynamic scaling of a boson Green's function). 
So such compressible states must have 



e = rf- 1. 



(1.6) 



Eq. (1.6) is a key constraint on holographic theories of compressible states with hidden Fermi 



surfaces, and our main results will be restricted to systems with IR metrics which obey (1.1) 



and (1.6). 



The arguments above also show that we can interpret Q as the dimension of the 
momentum-space surface on which there are zero-energy excitations. The value (1.6) is 
generic for systems with Fermi surfaces. 

We will compute the holographic entanglement entropy [55j, Se^ for theories obeying 



(1.1) and (1.6) for arbitrary smooth shapes of the entangling region(s), and general T in 
Section III At T = 0, and for a single connected entangling region as in Fig. [T| we find a 
logarithmic violation of the area law for = — 1, and to leading log accuracy 



5^ = ^ Q('^-i)/'iE ln(Q(<^-i)/'^E), 



(1.7) 



in the regime where the argument of the logarithm is large, and so S'^; can be computed 
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FIG. 1: Geometry of holographic entanglement. The d spatial co-ordinates are Xi = (xj^Xd), 
with j = 1, and the emergent holographic direction is r. The entangling region of the 

compressible quantum state is shown shaded; its boundary is the surface described in co-ordinate 
patches by x^ — w{xj), j = 1, ... ,0? — 1, which has surface area S. This boundary is extended 
along the holographic direction into the surface written locally as x^ = VF(r, Xj); the holographic 
entanglement is proportional to the surface area, A, of this extended boundary. The area A 
has to be minimized while keeping the surface x^ = '^(^j) fixed. 



entirely from scales where (1.1) applies. Here 77 is a dimensionless numerical constant which 
depends upon the couplings of the holographic theory, but is independent of Q and of any 
property of the entangling region. The only dependence on the entangling region is via S, 
which is its {d — l)-dimensional surface area (in = 2, S is the perimeter of the entangling 
region — see Fig.[T]). In general, gapless theories have an entanglement entropy which depends 
upon the full geometry of the entangling region [ 551457] . and not just on the surface area 
S. However, systems with gapless fermionic excitations on a (rf — l)-dimensional spherical 
Fermi surface [4314451 [581462] have an entanglement entropy which has precisely the form of 



Eq. (1.7), depending only on S and no other characteristic of the smooth entangling region. 

The specific Q-dependence of (1.7) arises from the Q-dependence of the numerical pref- 
actors of ( |1.1D . The latter dependence is quite complicated, and depends on details of the 
ultraviolet (UV) physics. However, all of this UV dependence cancels out in the entangle- 
ment entropy, and only the universal IR Q-dependence shown in (1.7) remains; this key 
result is established in Section II Upon using the Luttinger relation, Q ^ k^^ between the 
charge density and the Fermi wavevector fci?, the dependence in (1.7) is also identical to 
that found in the Fermi surface computations 



Thus, remarkably, the Lut- 
tinger relation of condensed matter physics is connected to some of the central principles of 
the holographic theory: Gauss's law [36j, the area law for entanglement entropy [55j, and 
the universal Q dependence of the metric of the holographic space. 
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We also examine the nature of entanglement between disjoint regions (the "mutual in- 
formation"), and again find agreement with the Fermi surface results. 



All of the above features of (1.5) and (1.7) support the hidden Fermi surface interpretation 
of holographic theories with the metric (1.1) for 9 = d — 1. 



We will also extend the EMD theories by adding explicit fermionic degrees of freedom in 



Section |IV[ These fermions are gauge-neutral mesinos, and form Fermi surfaces which are 
directly visible in the holographic theory. In the language of Ref . [29] , we are moving from 
the NFL state to a FL* state with co-existing quark (hidden) and mesino (visible) Fermi 
surfaces. We find that the FL* state also has the holographic entanglement entropy ( |1.7 ), 



but with Q-> Q- Qr, 



The FL* state has Qq^ark = Q - Qn 



and so this result is 



consistent with the holographic entanglement entropy being a measure of the hidden quark 
Fermi surfaces. 

We also find interesting inequalities for the values of z and 9 for a wide class of models 
with holographic duals, generalizing analogous relations obtained earlier [131 El ESI 140] . 
The considerations in Section [TTl] show that the 'area law' of the entanglement entropy in 
the IR (modulo logarithmic corrections) requires 



9<d-l, 



and so we expect this inequality to apply to generic local quantum field theories on the 
boundary. Note that the UV always makes a contribution that obeys the area law. Rela- 
tivistic conformal field theories, with = 0, do obey (1.8). However, there are supersymmet- 



ric lattice models [63j which do have large degeneracies in their spectrum, and these could 



violate (1.8). In Section II we find that the dynamic critical exponent obeys the inequality 



z>l + 



9 
d' 



:i.9) 



This inequality appears to be tied to the existence of a holographic gravity dual. It is quite 



remarkable that both inequalities (1.8) and (1.9) are realized as equalities ^40j only for gauge 
theories of non- Fermi liquid states in = 2, which have the exponents ^ = 1, z = 3/2 in 
analyses to three loops [53l[6l]- This is the first quantitative connection between the field- 
theoretic and holographic approaches to non- Fermi liquids. These results also suggest that 
there is a strongly coupled non- Fermi liquid state in = 3 with z = 5/3, in contrast to 
the weak-coupling field theory analysis [65j which only has marginal corrections to Fermi 
liquid theory. We summarize our comparison between the field-theoretic and holographic 
approaches to non- Fermi liquids in Table |T| 

The outline of our paper is as follows. In Section ITT] we will recast the EMD theories 



[IIlIISllIllIiniESlEalES] in a general context as reahzations of the IR metric (jLTj). We 
will pay particular attention to the Q dependence of the solution. Various aspects of the 



entanglement entropy of such theories will be studied in Section III, We will establish (1.8) 
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Field theory 



Holography 



A gauge-dependent Fermi surface of 
overdamped gapless fermions. 



Fermi surface is hidden. 



Thermal entropy density S ^ T^/^ in 
all d for hyperscaling violation expo- 
nent 9 = — 1, and z the dynamic 
critical exponent. 



Thermal entropy density S ^ T^/^ in 
= 2, where z is the dynamic critical 
exponent. 



Logarithmic violation of area law of en- 
tanglement entropy, with prefactor pro- 
portional to the product of and 
the boundary area of the entangling 
region. 



Logarithmic violation of area law of en- 
tanglement entropy for 9 = d — 1^ with 
prefactor proportional to the product of 
Q(d-i)/d ^Yie boundary area of the 
entangling region. 



Three-loop analysis shows z = 3/2 in 
d = 2. 



Existence of gravity dual implies z > 
1 + 9/d] leads to z > 3/2 for ^ = - 1 
md = 2. 



Fermi surface encloses a volume pro- 
portional to Q, as demanded by the 
Luttinger relation. 



The value oi obtained from the en- 
tanglement entropy implies the Fermi 
surface encloses a volume proportional 
to Q, as demanded by the Luttinger 
relation. 



Gauge neutral 'mesinos' reduce the 
volume enclosed by Fermi surfaces of 
gauge-charged fermions to Q — Qmesino- 



Gauge neutral 'mesinos' reduce the vol- 
ume enclosed by hidden Fermi surfaces 

to Q — Qmesino- 



TABLE L Comparison of the field-theoretic and holographic approaches to non- Fermi liquids 
(NFL). The field-theoretic results are as described recently in Refs. 



as a general inequality imposed by the area law of entanglement. For the hidden Fermi 
surface case, ^ = — 1, we will compute the T = entanglement entropy for an arbitrary 
shape, its crossover to the thermal entropy for T > 0, and the "mutual information" char- 
acterizing the entanglement of disconnected regions; all of these results will be of the same 
form as expected for systems with Fermi surfaces, with the specific Q dependence displayed 



in Eq. (L7). Section IV will extend the EMD theories to include gauge-neutral fermions, or 



'mesinos', which have visible Fermi surfaces. We will present a fully self-consistent compu- 
tation, including the back-reaction of the mesinos on the metric, using the approximation of 
the 'electron star' theories p!5l - [T7l [39] . We will find that the mesinos reduce the holographic 
entanglement entropy in precisely the manner expected from the overall Luttinger relation 
on the Fermi surfaces of the hidden quarks and visible mesinos [29j. Finally some general 



discussion appears in Section V 
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II. EINSTEIN-MAXWELL-DILATON THEORY 



We begin with a discussion of the basic characteristics of the EMD theory of compressible 
quantum states [HI [25] in d spatial dimensions. As reviewed in Ref. [37], the globally 
conserved U(l) charge of the compressible state is holographically realized by a U(l) gauge 
field Afj^. In addition, the EMD theories also include a scalar field (the 'dilaton') which 
is dual to a relevant perturbation on the UV conformal field theory, and which allows access 
to a wider range of IR scaling behavior in compressible states. 

So we consider the holographic Lagrangian 

defined on a (rf + 2)-dimensional spacetime with Ricci scalar i?, with Maxwell fiux 
associated with ^4^, and a dilaton field $ with potential V{^) and coupling Z{^). 

We use co-ordinates (t, r, x^), where t is the time direction, r is the emergent holographic 
direction, and Xi {i = 1 . . . d) are the fiat spatial directions. We will examine solutions with 
metric 

ds^ = {-f{r)de + g{r)dr^ + ^\ , (2.2) 



only the temporal component of the gauge field non-zero 

At = —hir), (2.3) 

K 

and a dilaton field $(r) dependent only upon r. Under these conditions, we can work with 
action per unit spacetime volume of the boundary theory 

dr—^\/J{r)g{r)CEMD- (2.4) 
The Einstein equations for this Lagrangian are 



dL^^Jir) V ^f{r) Sg{r) 

_ f'{r) ^ g'jr) ^ 4 ^ 4^'{rf _ ^ 
rf{r) rg{r) d 
2K''r\g{r)f'^ SSemd 
L^^W) Sg{r) 



(2.5) 



2r/(r) 2/(r) 2^^^^ 2r^ ' 
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while the equation of motion of the dilaton field is 



+ 



1 



2f{r) 4/(r) 
Finally, the only non-zero Maxwell equation is Gauss' Law, which yields 



(2.6) 



k^SSemd ^± h'{r)Zmr)) 
Ld Sh{r) dr \rd,/f{r)g{r) 



0. 



(2.7) 



The integration constant in (2.7) is set by the charge density on the boundary [36], and so 
we have 



(2i 



J rd^f{r)g{r) 

The dependence of the solutions on the charge density Q will be crucial to our purposes. 

We now discuss the structure of the solutions in the IR limit, r ^ oo. As we discussed in 
Section [l| we are interested in solutions which obey (1.1) in this limit. Extending the results 
of Ref. [25j to general rf, we can deduce that such a solution will emerge from the equations 
of motion provided we choose the large $ behavior to obey 



Z($) = Zoexp(a$) 

= -yoexp(-/3$) 



as $ ^ oc. 



(2.9) 



with a, P > 0, and the exponents z and 9 are then given by 
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a + {d- l)/3 

9_ 8{d{d-9) + 9f 
~^d~^ d^{d-9)a^ ' 



(2.10) 
(2.11) 



The inequality (1.9) is clearly obeyed by ( |2.11 ); this inequality can also be obtained by 
applying the null energy condition discussed by Ogawa et al [4Qj. Also, imposing the 



inequality (1.8) on (2.10), we obtain 



fd-i^ 



(2.12) 



{2d - 1) 

The compressible state with hidden Fermi surfaces has 9 = d — 1, and this requires that 



(2.12) is realized as an equality. 



Inserting (2.9) into the equations of motion, we can scale out the explicit dependence on 
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r, Q, Vo and Zq, by parametrizing the solution in the foUowing form 



fir) 
9{r) 
h{r) 



1 -2-2d{z-l)/{d-e) 



Q2/d Ql/d 



-2+2e/{d-e) 



V^-\VoZo)-'h 
V^-\VoZo)-'go 



Ql/d 



n -d-dz/{d-e) 



2d(l+0)/a 



where 



and 



e = 



Jd^ 

e 



d{d — 



(2.13) 



(2.14) 



(2.15) 



The metric (2.13) is of the form (1.1), as expected. 



To fully specify the IR solution, we need fix the values of the prefactors /o, go, ho and 



To in (2.13). An analysis of the equations of motion shows that the large r limit determines 
the values of 3 of these constants, fo/h^, go, and ro; the results take a remarkably simple 
form when expressed in terms of z and 9: 



= (z- 1)-' {z + d-e- if+^ (z + d-e) 



d^ 



fo ^ 

hi 



{d-ey 

{z - 1)-'-^ {z + d-e- 1)^+^ {z + d-e) 

{z-l){z + d-e -i)-^ 



(2.16) 



The solution for the hidden Fermi surface state is realized simply by setting ^ = — 1 in 
these general expressions. 



A key feature of our results in (2.13) and (2.16) is that all of the exponents, and some of 



the scale factors, are fully determined by the IR limit. In particular, and rg are fixed by 



(2.16) and so are independent of the charge Q. On the other hand, /o and are not fixed 
by the IR theory, only the ratio /o/Z^o fixed; these parameters will be set by matching 



to appropriate constraints arising from the underlying UV theory: see Section |IV A| and 
Ref. In general, the results of this matching can depend upon the value of Q, and so 
there can be additional Q dependence in our results via the values of /o and /iq, beyond that 



explicitly displayed in Eq. (2.13). However, this UV-induced Q-dependence does not infect 



the values of , '^o, and fo/h^^ and this will be crucial below in establishing the result (1.7) 
for the entanglement entropy. 

Gauss' law plays an especially important role in this result. Indeed, our arguments 
indicate that the bulk dual of the Luttinger constraint is eflFectively Gauss' law [36l|39j. As 
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far as the solutions are concerned, Gauss' law communicates a certain special combination 
of UV data to the IR in a way that depends only on the intervening visible bulk charge. 
The result is that the entanglement structure in the IR is sensitive to the UV only through 
the total charge density less the visible fermion charge, as we will see in more detail in 
Section IVl 



A. Non-zero temperatures 



The equations of motion (2.5 2.7) actually have a broader class of black hole solutions. 



associated with raising the compressible state to non-zero temperatures. These solutions 



modify the functions /(r) and g{r) in (2.13) to 

Mr) = f{r) (l - {r/ruf '^^/^'-''^^) 
9T{r) = g{r) (l - (r/r,f'^^^^'-'^^ 



(2.17) 



where Th is the position of the black hole horizon; the functions $(r) and h{r) remain 
unchanged. 

As usual, the requirement of the absence of a conical singularity at the horizon in Eu- 
clidean spacetime fixes the temperature T. To determine T we write r = Vh — p^^ and 
compute the metric in the limit p ^ 0; it has the schematic form 



Ci p^dt^ + C2 rfp^ + 



(2.18) 



where Ci^2 are p-independent constants. We now write t = i(j)/{27rT) and demand that the 
near-horizon metric have the planar contribution ^ p^dcf)^ + rfp^, which ensures the absence 
of a conical singularity provided (j) is periodic with period 27r. This fixes the Hawking 
temperature T = [c\jc^^l^ Such a computation yields the following relationship 
between r/^ and T: 



Tr 



dz/{d-e) 



Q 



-z/{d-e) 



ho{z + d-e) 



^ ^ 47r(z - 1) 

The thermal entropy density, aS, is given by the area of the horizon: 



(2.19) 



S 



2tv L'^ 

Qrp{d-e)/: 



2neL 



K 



4:Tr{z - 1) 
ho(z + d — 6 



(d-e)/z 



(2.20) 



and this is of the form in (1.4). Recall that the value of ho, and hence the prefactor in S, 



is not fully determined by the present IR solution. We have to embed the solution in a UV 
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AdS geometry to fix /iq, and this will be considered later in Section [IV| Such a determination 
can lead to a value of ho which is Q-dependent, as we noted earlier. So the Q-dependence 
of the thermal entropy is not only due to the factor of Q displayed explicitly in (2.20). This 
is compatible with the Fermi surface interpretation, because the dispersion of low energy 
excitations near the Fermi surface can have a complex Q dependence, and this does infiuence 
in the thermal entropy. In contrast, as will see in Section [TTTl the entanglement entropy has 
only the explicit Q-dependence shown in (1.7). 



III. SHAPE DEPENDENCE OF THE ENTANGLEMENT ENTROPY 

This section will extend the computation of the holographic entanglement entropy by 
Ogawa et al [40j to a general shape of the entangling region, to general and determine 
the full dependence on Q. We will find that the results correspond to those expected from 
a Fermi surface. 

For this computation, we will only need the equal-time metric, which we write as 

ds^ = ^{gy^/(d-o)dr^ + dx',) (3.1) 
where we initially allow for a general ^, and 

9o = Q''Vo-\VoZo)-'go. (3.2) 
We emphasize that the only Q dependence of go is that explicitly displayed above, because 



go is fully determined by (2.16). 

We consider an entangling region of an arbitrary shape. However, all our considerations 
will be associated with local area elements, and so we can always represent the boundary of 
the entangling region by solving for one of the spatial co-ordinates, say in terms of the 
remaining co-ordinates. So the surface is (see Fig. [T]) 

Xd = w{xj) , j = l,...,d-l (3.3) 

for some function w oi d — 1 co-ordinates. We will implicitly assume below that the indices 
j,/ = 1, . . . , — 1 below. The characteristic size of the entangling region, i?, is assumed to 
be much larger than the scale set by the density Q, and so kp^] under these conditions. 



the entanglement entropy can be computed using (3.1), and we don't have to include the 
crossover to the ultimate UV behavior near the boundary. 

The prescription for the holographic entanglement entropy is [55] to compute the minimal 



area of a surface which encloses this entangling region in the extended holographic space. 
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Let us parameterize the extended surface by (see Fig. [T]) 



Xd = W{r, Xj). 



(3.4) 



Then we have to find the optimum function W{r^Xj) subject to the constraint 



W{0, Xj) w{xj). 



(3.5) 



Let us compute the area of the general holographic surface in (3.4). The induced metric 
on this surface is 



da' 



L' 



dW 
dr 



, 2 ^dWdW , , dW dW 

dr + 2—-——drdxj + djj> + -— - — 
or oxj \ oxj oxji 



dx j dx j I 



(3.6) 

The area element on the surface is determined by the square-root of the determinant of the 
induced metric, which is 



— d'^-^x- 

^d-0/(d-e) " ^1 



1+ ^ r-'O/{d-0) (Q^^ 



dXa 



90 



dr 



1/2 



(3.7) 



We now observe that for d — 9/ {d — 9) > 1, which is equivalent to (1.8), the r integral 
is divergent as r 0: then the leading term to the integral over dA is an ultraviolet 
contribution proportional to S (see Fig. [T]) which yields the 'area law' of entanglement 



entropy. Thus we expect that the inequality (L8) applies to holographic duals of all generic 



local quantum field theories which do not have large accidental degeneracies in their low 
energy spectrum. Also, as we noted earlier, relativistic conformal field theories have 9 = 0. 

The remainder of this section limits consideration to the case 9 = d—1 oi interest in this 
paper, where we have a logarithmic violation of the area law. Let us study the nature of the 
r ^ limit more carefully. Let us expand W in this limit as 



W{r, Xj) = w{xj) + r'^a{xj) + . . . , r ^ 0, 



(3.^ 



where it remains to determine the exponent, n, of the leading correction, and a is an arbitrary 



function of the d — 1 co-ordinates. Inserting this in (3.7) we have 



dA = L''gl/'-d''-'x, 



dw\^ ^ „dw da r^^" '^^ n n 

I +2r" — — + — ^ nV + .. 



dxjj 



dxj dxj 



9o 



1/2 



(3.9) 



The variational derivative of the integral of this expression with respect to o'{xj) must 
vanish. A non-trivial solution is only possible if the two leading terms in powers of r can 
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cancel against each other. So we must have n = 2{n — d) or 



n = 2d. 



(3.10) 



So the r- and a-dependent terms inside the square-root in are indeed subdominant, and to 
leading logarithmic accuracy we can write 



Se — 



27r 



dA = 



1/2. 



dr 
r 



where 



E = d'^ 



1 



dw 

dXj 



-,1/2 



(3.11) 



(3.12) 



The quantity S depends only on the entanghng region on the boundary, and indeed it is just 
its surface area. So we conclude that the log-divergent entanglement entropy is proportional 
to the surface area of the entangling region, and is otherwise independent of its shape. This 
is precisely the property of the entanglement entropy of a spherical Fermi surface [431 [58] : 
our holographic analysis is for spatially isotropic systems, so a spherical Fermi surface is 



expected. Also note from (3.2) that the prefactor of (3.11) is of order Q^^ and so the 



complete Q-dependence of the entanglement entropy is that displayed in < \1.7\ i. 

Ogawa et al. [40] presented computations of the entanglement entropy in = 2 for two 
choices of the entangling region: a strip and a disk. Their result for the strip agrees with our 



general result (3.11), and also with the exponent value in (3.10). They presented numerical 



results for the disk, and the prefactor of their logarithm equals that predicted by (3.11). 



We also wish to point out that this result may be interpreted as additional evidence 
for the Ryu-Takayanagi formula. While their formula has been proven for spheres [67j, it 
remains unproven in general and is known to be modified in higher derivative gravity. In 
general, we can only show that the Ryu-Takayanagi proposal has the right basic structure to 
give an entanglement entropy. Our calculation here shows that the Ryu-Takayanagi formula 
reproduces in detail a universal feature of the entanglement structure of compressible states 
as expected from field theory for all region shapes. 



Finally, let us discuss the limits on the r integration in (3.11). We expect the large r 



limit to be set by the size of the entangling region, R. From ( |1.2[ ) and (1.3), we see that for 
9 = d — 1^ scales as Xj ^ R. So we can expect kFrmd.^ ^ {kpRY^^. We will discuss the 
value of Tmin more carefully in Section [IV| where we will argue that kpr^^^ ^ 1. Note that 
^max ^ ^min bccausc kpR ^ 1. With these limits, and using Q ^ kp and S ^ R^~^^ 



we see 



that (3.11) is of the form (1.7) 
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A. Entanglement to thermal crossover 



Using our finite temperature solution, we may also study the von Neumann entropy of a 
subsystem at finite temperature. There will be a universal crossover function that connects 
the anomalous entanglement entropy to the thermal entropy discussed above. This crossover 
function is universal because the same low energy degrees of freedom responsible for the 
long range entanglement also control the low temperature thermodynamics. However, this 
universality is only defined up to boundary law terms since these may be generated by UV 
processes and can depend on the UV cutoff. On general scaling grounds, this crossover 
function will depend on region size R and temperature T according to 

S{R, T) = Q^^-^)'^T^Fs[RT^''). (3.13) 

Note that we have suppressed a dimensionful constant multiplying RT^^^ from the dispersion 
relation analogous to the Fermi velocity of a Fermi liquid. The scaling function itself is 
independent of this constant. A general discussion of such scaling forms will appear very 
soon in [68j where it is shown that = (1 — for a compact region. This result is obtained 
by matching onto the thermodynamic entropy where extensivity requires Fs{x) x^. We 
may obtain this function holographically by computing areas of bulk minimal surfaces in 
black hole backgrounds. We set = — 1 in this section. 

We consider a strip-like entangling region of cross section S/2 and width R. Assuming 
jld-i ^ ^ then the total surface area is S. Because of the translation invariance in the 
cross-sectional directions, the bulk coordinate r will be a function of only one boundary 
coordinate r = r{x). The area is 




^ = (?) dx(-\ Jl + goT^i<^-^){l-{TlTHf'^^A^^ (3.14) 



-R/2 



dx J 



and the resulting equation of motion has an integral of the motion given by 

1 1 



1 

max 



Solving for dr /dx we find 



(3.15) 



dr V^ma.^-r^'^ 



(3.16) 
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We now recast the area as an integral over r to obtain 



and 



R = 2 dr =■ (3.18) 

We may safely set the UV cutoff Tmin = in the integral for R. This gives a universal 
relationship between i?, Tmax, and r/^. Changing variables to u = r/r^g,^ we find 

R = 2^orLl du ; (3.19) 



According to (2.22) we have ^ T~-^/^ and upon putting 2v^r^ax = RG{RT^^^) we find 
1 = G{RT^'') / rf^ , = (3.20) 

where we have absorbed a dimensionful constant into RT^^^ . This constant can depend 



on Q as well as other UV data (through /iq as in (2.19)) just as the Fermi velocity in a 
Fermi liquid can, but as we have shown explicitly here, the form of the scaling function is 
independent of this constant. 

We now determine the function G in certain limits. Starting from 

1 = Gix) / dii. (3.21) 
Jo Vl - ^'Vl - ud^^+^){xG{x)Y+^ 

we see immediately that G{x) goes to a constant as x ^ 0. On the other hand, as x ^ oc 
we have G{x) = 1/x to keep the integral well defined. This immediately gives r^^^ ^ 
as expected. To compute the subleading term at large x we define G{x) = {1 — 5)/x which 
gives 

1 = - / du ^ , (3.22) 

xj^ Vl - u^^^Jl - 2i^(i+^)(l -{z + 1)5) 

There is a logarithmic singularity as 5 ^ coming from 2i ^ 1, so to extract it we set 
u = 1 — V and expand 

rO{i) I 
X = dv , — (3.23) 

Jo V2dv^d{l + z)v + {1 + z)S 
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This gives 

x^2d^{l + z) - In (l/S) (3.24) 

and hence 5 ^ exp {-d^2{l + z)x) with X = RT^'\ 

Having estabhshed a scahng form for Tmax, a similar form immediately follows for the 
area itself. We have 

A=(^\ 2y^L^ f ^ ^ ^ ^ (3.25) 

Using a/^ ~ 

we see that the area reduces to 

Q(rf-i)/rfSln(Q^/^i?) (3.26) 

for RT^'^ and to 

Q{d-i)/d^-^j.i/z ^3^27) 

for RT^/' oc. 

The crossover function has = instead of (j) = {1 — d) / z because we worked in a limit 
where the cross-section S was larger than any other scale. In that limit matching to the 
thermodynamic result demands Fs{x) x instead of Fs{x) x^ and hence = 0. 

The interpretation of this result in the field theory is simple: we have a rf — 1 dimensional 
manifold (the Fermi surface) of gapless modes propagating in one eflFective (radial) direction 
with dynamic exponent z. At finite temperature these modes give an entropy of the form 

zero temperature each mode gives a In (Q^^^R) contribution to the 
entanglement entropy for a total entropy of Q(^~^)/^i;in (Q-^/^i?). 



B. Mutual information 

We can also compute the entanglement entropy for disjoint regions using the Ryu- 
Takayanagi formula [55]. Such configurations are interesting in part because they can be 
used to define subtracted versions of the entanglement entropy that are not UV sensitive. 
For two distant regions we find that S{AB) = S{A) + S{B). As the regions are brought 
close, the bulk minimal surface experiences a phase transition and suddenly connects the 
two regions in the bulk. Using S{A)^ S{B)^ and S{AB) we may compute the mutual infor- 
mation /(A, B) = S{A) + S{B) — S{AB). This quantity is guaranteed to be positive and 
has a variety of other favorable properties. It is also insenstive to the UV and measures 
both quantum and classical correlations. 

For two strips of cross-section S/2, width and separation X, and with Q~^/^ <C X <C 
W (see Fig. [2]), we computed the holographic entanglement entropy from (1.1) for ^ = rf — 1, 
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-XI2-W -X/2 X/2 XI2 + W 



FIG. 2: Geometry of mutual entanglement bewteen regions A and B. 



and found a mutual information of the form 



(3.28) 



which is identical to the form obtained from a Fermi surface [69]. For some critical value 
of X of order W ^ the mutual information will drop continuously to zero and remain there 
as X is increased. In principle, the long distance decay of the mutual information should 
contain information about the decay of correlation functions, but the geometric prescription 
is blind to this subleading contribution. 



The geometric computation proceeds from (3.17) and (3.18) which give the area for a 



minimal surface of boundary width R. We set r/^ = oc in this section. Region A is a strip 
between x = X/2 and x = X/2 + W and region 5 is a strip between x = —X/2 — W and 



X = —X/2. The entropies S{A) and S{B) are identical and proportional to (3.17) with 
R = W. For X <C the entropy S{AB) consists of two pieces, one associated with a 
surface connecting x = —X/2 — Vl^ to x = X/2 + W and the other associated with a surface 
connecting x = —X/2 to x = X/2 (see Fig. [2]). These surfaces both obey (3.17) and (3.18) 
but with R = X + 2W and R = X respectively. 

Setting r/^ = oc in (3.17) and (3.18) and changing variables to r = r^ax'^ we have 



A={-] (2LV^o) 



du 



U 



2d 



and 



R = {2^orL.) [ 



2d 



du u 



u 



2d 



(3.29) 



(3.30) 
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Taking = e ^ we find 

^= (f)(2L''V^)ln(^) (3.31) 

and 

R = 2^,rij, (3.32) 
witli k = duu^^~^/Vl — v?^. Combining these two equations immediately gives 

^W^(f)(2L^VS)ll„(^). (3.33) 

The mutual information is 

/ (X A{W) + A{W) - {A{X) + A{X + 2W)) (3.34) 

which gives 

(|),2i^^)l,„(_Jf!_). (3.35) 

Note that the UV cutoff e has vanished and that in the limit X <^W the argument of the 



log simplifies to W/X as claimed in (3.28). 



IV. ADDING MESINOS 

All the analysis of this section will be restricted to spatial dimension = 2, and 9 = d—1. 

This section will add fermionic matter to the NFL phase of the EMD theory, in a manner 
analogous to previous computations [21 HHHl [MEl [IMll [20l[25l|30l[3ll[35l[36l[^ 
These fermions represent gauge neutral 'mesinos' on the boundary, in the language of 
Refs. [36i [39j . When a non-zero density of such mesinos is present, we are in a FL* phase, 
and there is a visible Fermi surface, which encloses a volume associated with the mesino den- 
sity, Qmesino- Wc are primarily interested here in the influence of this visible Fermi surface 
on the hidden Fermi surfaces. 

The mesinos occupy fermionic states obeying the Dirac equation which can be written in 
terms of a two-component spinor, Xk [SI ESJ [36j : 



-^^Ty^|:-mLa-y7H-fc^Vy7(;0-^A(Ojx/c(r) = ^X/c(O (4-1) 

where q is the conserved Q-charge of the mesinos, a are Pauli matrices acting on the spinor 
space, k is the momentum of the state along the boundary, and E is the energy eigenvalue. 
The mass m determines the scaling dimension of the mesonic operator on the UV conformal 
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field theory, and we will impose the UV boundary condition ^ 0) ^ r^^, and outgoing 
wave restrictions for r ^ oc. The action associated with these fermionic states can be 
obtained by evaluating the determinant of the Dirac operator, as discussed in Ref. [7J, which 
is expressed as a sum over quasi-normal modes. This is a rather involved computation, 
but we can obtain some features of interest to us by using the structure of the solutions 



of (4.1), which were discussed in some detail by lizuka et al. [25j. For the metric (1.1), 
we are in the "Fermi liquid" regime examined by lizuka et al: the quasi-bound states are 
well-defined near the Fermi energy, ^ 0, with an exponentially small decay rate which 
can be computed by WKB methods. For now, let us ignore the tunneling decay process. 
Then we have discrete set of energy eigenstates labeled by the discrete index ^, just 

as in the simpler situation discussed in Ref. [36j. The action of these fermionic states is 
obtained by adding up the energy of all the occupied negative energy states [71 [36] 

.S^esmo = E / S ^^(^) ^ i-Ei{k)) (4.2) 

If we included the decay processes, this expression would be replaced by the determinant of 
the Dirac operator [7J . 

We are now in a position to write down the corrections to the equations of motion 



(2.5), (2.6), (2.7) due to the presence of the mesinos. The equations now follow from the 



stationarity conditions on the complete action Semd + ^mesino- The functional derivatives 



of Semd were given in (2.5), (2.6), (2.7), while those of 5mesino can be evaluated using the 



Feynman-Hellman theorem ^j: we normalize the fermion states so that 
and then the functional derivatives of (Smesmo follow from those of Ei{k): 



- xUr) C^^i) (4.4) 



Sgir) ' V 2 g{r)dr. 

For our purposes, the most important modification is to Gauss' Law, which is modified from 
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(2.8) to 



L h\r)Z{^{r)) 
Kc r^^f{r)g{r) 



4^ 




e. (4.5) 



These equations are too difficult to solve in general, but we can easily identify the char- 
acteristic scale of r at which the mesino states are present. And this follows by generalizing 



the dimensional arguments leading to (2.13). Using the Q dependence shown explicitly in 



(|2.13|) , we deduce that the Q-dependence of the Dirac wavefunction can be written in the 

Xk{r) = Fk{QT^) (4.6) 



scaling form 



where we reiterate that we are in = 2. Also, recall that below (2.13) we discussed additional 



Q-dependence, arising from the UV matching, in the values of /o and /iq; examination of the 
structure of (4.1), while using the fact that /o/Z^o Q-independent, leads to the conclusion 



that this additional Q-dependence does not infect the scaling form in (4.6). It does, however. 



appear in the value of the eigenvalue E. From these results we conclude that all the occupied 
mesino states are at r ^ ^/Q. 

A similar conclusion is reached by a WKB analysis of the Dirac equation, along the lines 
of Ref. [17j. By an appropriate change of variables, the Dirac equation is converted to a 
Schrodinger equation for the top component of the wavefunction, and the mesino wavefunc- 
tion is then concentrated in the classically allowed region of the potential: this argument 
also leads to a characteristic r ^ \fQ. 

We now look at the regime r ^ \fQ which is relevant for the computation of the en- 
tanglement entropy. The mesino states are at a scale r ^ a/Q, and we will ignore their 
possible contribution from tunneling out of the classically allowed region into r ^ \fQ. The 
possibility remains that these tunneling processes have a strong back-reaction on the metric, 
but we will not study this here. With this assumption, we examine Gauss' Law for r ^ \fQ\ 



we can write (4.5) as [391 



L /iXr)Z($(r)) 



r > VQ. 



(4.7) 



where, after using (4.3), we have 

Qmet 



47^2 



9{-Ee{k)) 



(4i 



is the total charge density in the occupied mesino states. So for r y Q, (4.7) replaces (2.8), 



while the differential Einstein equations (2.5) remain unchanged. Also, there is no mesino 



contribution to the dilaton equation (2.6), which remains unchanged. So the solutions for 
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r ^ ^/Q are as in (2.13) and (2.16), but with Q ^ Q — Qmesmo in the exphcitly displayed 
Q dependence. This is the key result of this section so far. 

Also note that (2.5), (2.6) and (4.7) do not determine the values of /o and ho. So, in 
general, the presence of the mesinos will change the values of /o and /iq, beyond that associ- 
ated with Q ^ Q — Qmesino- Howcvcr, go was fully determined by (2.5), (2.6) and (4.7), and 
this is the only parameter that appears in the holographic entanglement entropy. So we can 
simply apply the computation of the holographic entanglement entropy of Section III after 
the substitution Q ^ Q — Qmesmo- The present argument also shows that this computation 
applies for r ^ ^/Q, and so Tmin ^ a/Q- We schematically illustrate this result in Fig. [3| 

So, as claimed earlier, we have established that the main effect of the occupied mesino 
states in the FL* phase is to subtract away from the charge appearing in the holographic 
entanglement entropy. This is consistent with the expected transfer of quarks to mesino 
states, and the corresponding expected decrease in the hidden Fermi surface volume of 
the quark states. This feature provides additional support to the hidden Fermi surface 
interpretation of the NFL state. 



A. Thomas-Fermi theory 

Our analysis of the Dirac equation above treats the mesino states in a Hartree-Fock-like 
theory. The resulting equations are quite complicated to solve, and so we will consider here 
the analog of the simpler Thomas- Fermi theory. In the present context, the Thomas- Fermi 
theory is the 'electron star' approach p^SHTTl [39] . in which the electrons are treated as a 
continuous fluid which obeys the equation of state of free Dirac fermions in a local chemical 
potential. This fluid resides in the classically allowed region of the potential associated with 
the Dirac equation fTT]. In this approach, we will be able to embed our IR solution in an 
EMD theory which has a AdS4 metric in the UV. We will compute the transition from the 
NFL state to a FL* state with occupied mesino states, while including the back-reaction 
of the mesinos on the metric. We will also flnd here that the mesino states appear at a 
characteristic r ^ ^/Q. 

In [39] is was shown that one can study charge fractionalization in the EMD theories 
by dialing a relevant coupling in the dual UV theory. This coupling is dual to the dilaton 
field and can be made relevant by choosing the dilaton mass appropriately. The analysis 



presented in this subsection will closely follow [39j. We consider the EMD theory (2.1), and 
add the contribution from the Dirac fluid, which is simply given by the pressure 

Auid =p{l^\oc)' (4.9) 
Here /xioc stands for the local chemical potential, which for irrotational, zero temperature 
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FIG. 3: Schematic illustration of the holographic geometries of various compressible phases. The 
phases are labeled non- Fermi liquid (NFL), fractionalized Fermi liquid (FL*), and Fermi liquid 
(FL), following the notation of Ref. [29]. The boundary theory at r = has total charge density 
Q which sources the bulk electric field £r] we define £r to equal the left-hand-side of (2.8). The 
NFL phase is described by the theory in Section [llj and has the IR metric (1.1). The shading 



in the bulk region represents the density of mesinos. The FL* phase is described in Section IV A 
using the Thomas-Fermi approximation. The FL* phase has mesinos are in the bulk shaded region 



at r ^ vQ, which this corresponds to the shaded region in Fig. [5j the region r ^ \/ with 
£r — Q — Qmesino, determines the entanglement entropy. The FL phase is not described in the 
present paper: its geometry is confining and terminates at a finite r where £r — 0, as discussed in 
Ref. 
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fluids can be taken to be 



e h 



(4.10) 



where we used (2.2) and (2.3). The pressure of the fluid can be expressed in terms of the 



h 



energy and charge density 

— p = p T=(T, a = f3 eV — m?de^ p = f3 j e^V — rh?de^ 

V J Jrh 

where we introduced dimensionless variables 



(4.11) 



1 



P 



P: 



G 



eL'^K 



a . 



K 



-m 



(4.12) 



The equations of motion for the Einstein-Maxwell-dilaton-charged fluid theory for the metric 



and vector potential ansatzes (2.2) and (2.3) are found to be 

1 g' 4\ gh ^ 



0, 
0, 



2 V / % + 4 4/ 



(4.13) 



One easily verifles that the above equations reduce to the EMD equations of motion (see 



Section II) in the absence of the fluid, that is when /iioc < rh. As in the previous section 
we see that the presence of the fermions in the bulk modifles Gauss' Law. In the fluid 
approximation the total charge carried by the fermions is 



(4.14) 



n 



where r\ and r2 > r\ are the radii between which a non-zero density of the fluid is present. 
The total charge now reads Q = Qmesino + Qquark, where 



quark 



Zmr))h'{r) 



for r > r2- 



(4.15) 



To explore the Q dependence of theory, it is useful to note here that the equations of motion, 
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(4.13), are invariant under the scale transformations 



r 




Xr 


/ 




f 


9 




X-'9 


h 




h 








2quark 







where A is the rescahng factor, and the additional rule that {5-,p,j5} {^,P,j5}. 

To embed the theory in an asymptotic AdS4 metric, we choose the dilaton potential to 
take the form near $ = of 



= -6 + 2M|L^$^ + (4.17) 

For computational convenience we will take the dilaton mass to satisfy M| = — such 
that the dual operator O has scaling dimension A = 2. Note that the mass is above the 
Breitenlohner-Freedman bound. If we also choose Z{0) = 1, the near boundary expansion 
of the dilaton field takes the form 

$ = (/)or + ^r^ + ... (4.18) 

where 00 and {O) are the source and the expectation value of the relevant operator O dual 
to the bulk dilaton $. 

It is now convenient to fix V{^) and Z{^) such that they interpolate between the AdS4 
solution in the UV and the IR solution described in Section HIl Note that indeed for the 



latter solution, (2.13), the local chemical potential goes to zero in the IR, which implies that 
the fiuid is absent. We expect our conclusions to be insensitive to the specific interpolation 
form, and we found it convenient to use the following expressions 

V{^) = -t; ^r-^n^ + (1^0 - eVl - tanh(a$/3)'). 



2cosh(a$/3) V2 
Z($) = exp(a$), (4.19) 

with Vo = 24 (a^ + 6) /a'^ to get the desired dilaton mass in the UV. One readily verifies 



that these expressions take the form of (2.9) with = 2 in the limit of $ ^ oc in the IR. 
Finally, without loss of generality we set a = 3 for computational convenience. 

For the full potential, the IR solution is modified by a series expansion together with 
a perturbation that will allow the fiow to be integrated up to the UV. Furthermore, the 



generalized Gauss' Law implies that we have to replace Q by Qquark in (2.13). The IR 
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solution thus reads: 




^quark 

9 = |e,„„.fi + E|; + *f"l' ("-21) 

\ n=l 



^quark \ n=l 



$ = 21og(V^r)+log(^iy|j +f:^ + '^^-'^- (4-23) 

The coefficients {fm 9n^ hn^Pn} are uniquely determined by the equations of motion, while 
the perturbation {5/, dg^ dh^ dc/)} has an overall free magnitude. This irrelevant perturbation 
has 

A^ = -- (V79-3) < 0. (4.24) 



2 

Different magnitudes for the perturbation will map onto different values of the dimensionless 
ratio of relevant couplings 00 //i in the UV theory, where /i is the boundary chemical potential 
rescaled by appropriate factors of e and L. 



We solve the equations of motion, (4.13), by numerically integrating up the IR solution. 



(4.23), to the UV boundary at small values of r. The only free parameter is the strength 



of the irrelevant perturbation in the IR. We can read off the relevant coupling 00 near the 



boundary by fitting to the boundary expansion of $ (4.18). Similarly, the boundary chemical 
potential, /i, and the total charge, Q, can be found from fitting to the boundary expansion of 
h{r) = c{jl — Qr) + . . . . A key feature of this IR-UV matching procedure is that the several 
parameters characterizing the UV limit to do not modify the IR limit of and hence do 
not infiuence the Q dependence of the entanglement entropy, as we have also emphasized in 
Section |l| In performing the integration we must set {a^p^p} to zero whenever /iioc(^) < ^ 
and keep them in the equations otherwise. 

In figure [4| we plot the ratio of the charge carried by the quarks and the total charge as 
a function of the relevant coupling, 00. As in [39j there is a third order transition between 
a partially fractionalized phase and a fully fractionalized phase. In the fully fractionalized 
phase we have ft\oc{r) < rh throughout the entire bulk. It follows that the total charge 
Q = Qquark- the partially fractionalized phase the fiuid is present in a region between ri 
and r2 in the bulk. In this phase both quarks and mesinos contribute to the total charge: 

Q = Qquark + Q mesmo • 

We have argued throughout this paper that for the present Einstein-Maxwell-dilaton-fiuid 
theory the holographic entanglement entropy indicates the presence of hidden Fermi surfaces 
of gauge-charged quarks in both phases. Furthermore, the presence of Fermi surfaces of 
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FIG. 4: Ratio of fractionalized charge to total charge Qquark/ Q as a function of the relevant coupling 
00 for {m, /3} = {3/16, 10}. The black dot denotes the location of the transition between the FL* 
and NFL phases. In the FL* phase indicated by the blue drawn line Qquark/ Q < 1, whereas in the 
NFL phase indicated by the dashed red line we have Qquark/ Q = 1- 



gauge-neutral mesinos associated with the fluid was established in [T6] . We can thus identify 
this transition with a phase transition between FL* and NFL. It is clear from (4.23) that the 
entanglement entropy in the FL* phase is also given by (1.7), but with Q Qquark- Via the 
Luttinger relation we flnd that this is consistent with the interpretation of the holographic 
entanglement entropy being a measure of the hidden quark Fermi surfaces. 

In flgure[5] we show the region in the bulk where the fluid is present; see also Fig. |3) It 
is clear that this region shrinks to zero at the FL* to NFL transition. It follows from the 
scaling transformations (4.16) that /iioc(^) = ^(Qquark^^)- From this we conclude that the 
fluid is present at r ^ y^Q^^mrk- The discussion below (4.6), applied to the present context, 
implies that the UV matching conditions do not modify this conclusion. 

For completeness we mention that as in [39] the present theory has two more IR solutions. 
These solutions describe the phase diagram beyond the FL* phase. There is a Lifshitz 
solution where the dilaton is constant. This solution can be reached from the UV by flne 
tuning the coupling to the relevant operator, 00, and bounds the partially fractionalized 
phase. Finally, there is a solution with $ ^ — oc in the IR. In [39j this was a domain 
wall-like solution plus a fluid corresponding to a phase where the flux vanishes in the IR and 
all the charge is accounted for by mesino Fermi surfaces. For the present theory, however, 
the solution is quite different. In particular, the fluid density diverges in the IR and so does 
the flux. We leave the study of this phase for future work. We just mention that this phase 
is not present if one chooses Z symmetric under $ — $, Z{^) = cosh(a$) for instance. 
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FIG. 5: The plot shows the minimal and maximal radii for which the fluid is present as a function 
of the relevant coupling (po. The factors of ft ensure that all quantities are given in dimensionless 
units. Remember that the IR corresponds to r — > oo. The black dot denotes the location of the 
transition between the FL* and NFL phases. The fluid is present at all radii for which ftiodr) > m, 
this is indicated by the shaded region; this corresponds to the shaded region in Fig. |3| It follows 
from scaling arguments that the fluid is present at r ^ y^Qquark- 




V. DISCUSSION 



We conclude by discussing some issues raised by our results for future work. 

The holographic entanglement entropy formula of Ryu-Takayanagi appears to detect 
only the hidden Fermi surfaces. The visible Fermi surfaces of mesinos in our work, and in 
previous work [21 HHBl [HI EU |36l [39] , do properly reduce the net charge associated with 
the hidden Fermi surfaces, but do not contribute themselves to the entanglement associated 
with the Ryu-Takayanagi result. Presumably, their contributions will appear in fluctuation 
contributions to the entanglement entropy: it would be useful to sort this out in future work. 

We have analyzed zero temperature solutions of the EMD theory that possess logarithmic 
violations of the boundary law for entanglement entropy. We now make a few additional 
comments on the singularity structure of these solutions along the lines of Ref. [48j (see also 
Refs. [131 [13 )• We focus on the case of = 2 for concreteness, but our remarks apply more 
generally. As a preliminary, we note that the components of the curvature tensor are all 
flnite in the coordinate system of (2.2) for / ^ j^-^z-i) g ^ diS in (1.1). However, 
the spacetime is not geodesically complete. Consider a radial geodesic with tangent vector 
n = (t, r, 0, 0) where r indicates a derivative with respect to the proper time r. By imposing 
the normalization condition n-n = —1 and using the existence of a conserved "energy" n • dt 
we can flnd a simple diflFerential equation for r(r). Omitting the details, we flnd an equation 
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of the form 




1 



) 



(5.1) 



with To and tq some constants. For large r this equation may be integrated from to yield 



This geodesic reaches r = (X) in finite proper time Too given by 



and hence the spacetime is not geodesically complete for z > 1. The inequality (1.9) with 
9 = d — 1 guarantees that our solutions always have this finite proper time singularity. We 
also expect to develop tidal singularities analogous to those in [48] in a parallel transported 
frame. These observations have practical importance in that we should apply absorbing 
boundary conditions at r = (X) when computing correlation functions. 

We also have perfectly sensible finite temperature solutions that still display a logarith- 
mic violation of the boundary law for low enough temperatures. One possible interpretation 
of this situation is to regard the singular zero temperature solution as a kind of interme- 
diate scale description of the physics that, however, may not be extendible all the way to 
zero temperature. A small finite temperature provides a regulator and these solutions still 
describe physics in a wide range of energy scales. In Ref. [40j, various other spatial metrics 
were written (without full solutions) that give a greater than logarithmic violation of the 
boundary law. From the perspective of field theory, these solutions are quite anomalous in 
terms of their entanglement entropy, so like the extremal Reissner-Nordstrom black hole, we 
believe that such a description must not be valid to arbitrarily low energies. As argued in 
Ref. 08] , such a breakdown is visible on the gravity side in that the deep IR geometry tends 
to excite even very heavy bulk degrees of freedom. 

There is an important distinction between these more anomalous states and the hidden 
Fermi surfaces we have considered: the hidden Fermi surface system can be stable to zero 
temperature. In fact, this stability is a non-trivial question on the field theory side especially 
given the propensity of the Fermi liquid to form paired superfiuid states. For a Fermi surface 
coupled to a U(l) gauge field, there is a range of parameters where the system is stable to 
pairing [TOj. However, the situation for a Fermi surface coupled to a non-Abelian gauge field, 
as may be more appropriate in the EMD context, is less favorable. The colored "quarks" may 
be paired into a wider variety of color states owing to the non-Abelian nature of the group, 
and invariably at least one of these channels is attractive. Whether this conclusion can be 
modified by the inclusion of other strongly coupled IR degrees of freedom, such as those that 
appear in the EMD theory, is not known. Thus although there remain issues associated with 




(5.2) 
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the singularity structure of our solutions, the physics they appear to describe is compatible 
with a true zero temperature compressible state. In the future, it would be interesting to 
develop additional bulk criteria that can explicitly rule out the anomalous solutions. 

We mention one final interesting point. Free fermions and Fermi liquids are also known to 
have anomalous fiuctuations in conserved charges ^58^ |69j . Such fiuctuations are interesting 
in part because of their greater experimental accessibility. Let Qa = Q be the charge 

operator restricted to region A of linear size R. Even if the total charge Q commutes 
with the Hamiltonian, the restricted operator Qa need not. It can be shown that CFTs 
with a conserved charge have {Qa) = and fiuctuations AQ\ = {Q\) ^ R^~^ (rf > 1). 
The compressible Fermi systems mentioned above have instead {Qa) and /S.Q\ = 
{{Qa — {Qa)Y) ^ R^~^\og{R). We suspect similar charge fiuctuations will be present 
in general compressible states, so it would be quite interesting to perform a holographic 
calculation of these fiuctuations. 
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